Superconductivity from Undressing. II. Single Particle Green's Function and 

Photoemission in Cuprates 



J. E. Hirsch 

Department of Physics, University of California, San Diego 
La Jolla, CA 92093-0319 
(February 1, 2008) 

Experimental evidence indicates that the superconducting transition in high Tc cuprates is an 
'undressing' transition. Microscopic mechanisms giving rise to this physics were discussed in the 
first paper of this series. Here we discuss the calculation of the single particle Green's function and 
spectral function for Hamiltonians describing undressing transitions in the normal and supercon- 
ducting states. A single parameter, T, describes the strength of the undressing process and drives 
the transition to superconductivity. In the normal state, the spectral function evolves from predom- 
inantly incoherent to partly coherent as the hole concentration increases. In the superconducting 
state, the 'normal' Green's function acquires a contribution from the anomalous Green's function 
when T is non-zero; the resulting contribution to the spectral function is positive for hole extraction 
and negative for hole injection. It is proposed that these results explain the observation of sharp 
quasiparticle states in the superconducting state of cuprates along the (7r,0) direction and their 
absence along the (7r,7r) direction. 



I. INTRODUCTION 

Photoemission and optical experiments indicate that in high Tc cuprates a transition from an incoherent to a 
partially coherent state occurs both as the hole doping increases in the normal state and as the system goes su- 
perconducting 1^-^. Basov and coworkers [Q have observed a lowering of c-axis kinetic energy as the transition to 
the superconducting state occurs in several cuprates, especially in the underdoped situation. It has been established 
however that the magnitude of the c-axis kinetic energy lowering detected is far too small to account for the supercon- 
ducting condensation energy at least in some cuprates ||]. Ding, Campuzano, Norman and coworkers [QjlJ^ as well as 
Feng and coworkers j|] have reported observations of sharp quasiparticle peaks in the superconducting state in angle 
resolved photoemission emerging from a highly incoherent normal state background along the (tt, 0) direction, close to 
the {tt /a, 0) point. Ding and coworkers have interpreted the photoemission peak in terms of an enhanced quasiparticle 
weight Z in the superconducting state, and Feng and coworkersH] have suggested that the peak in photoemission 
is a signature of the superfluid density. Norman and coworkers [[7|| have analyzed the photoemission observations in 
terms of a 'mode model' and emphasized the close connection between their observations and Basov's observation of 
kinetic energy lowering. Furthermore, Basov and coworkers have emphasized that kinetic energy lowering seems to 
occur only when there is a high degree of incoherence in the normal state, and appears to vanish as the normal state 
becomes more coherent (overdoped regime) |^ . They have furthermore proposed that the photoemission experiments 
suggest that kinetic energy lowering may occur also in — plane in the cuprates albeit only along the (tt, 0) direction, 
and for that reason may be difficult to observe directly. 

The model of hole superconductivity |^,^ predicted, before the experimental observations, that the superconducting 
condensation energy originates in in-plane kinetic energy lowering [ pO[ and arises from a process of undressing of hole 
^ ' carriers as the pairing state develops [pd[| . Thus it describes both the kinetic energy lowering, arising from the low 
• I— I . energy effective Hamiltonian, as well as the high energy optical spectral weight transfer, that has also been observed 
■ experimentally . In the first paper of this series (hereafter referred to as I) wc formulated more generally the 
^ , principles of superconductivity through hole undressing, and pointed out that this physics would show up both in 
■ - - ' the one and two-particle Green's functions, in qualitative agreement with the observations reported above. Here we 
report calculation of the single particle Green's function and spectral weight in the superconducting state and discuss 
implications for the understanding of photoemission experiments. 



1 



II. GENERAL PRINCIPLES 



In the class of models discussed in I, the wave function renormalization of quasiparticles is a function of the site 
occupation in a local representation. The wavefunction renormalization arises from coupling to a local boson degree 
of freedom. Three examples of specific microscopic Hamiltonians describing this physics were discussed in I. The 
'coherent' part of the electron creation operator at site i is defined by the following transformation: 

dl = [T-{T- S)j\u.-a]dV (1) 

with < 5 < r < 1. The J operators in Eq. (1) are quasiparticle operators 10, and n^j is the electron site occupation. 
Eq. (1) expresses the fact that the electron becomes less coherent as more electrons are added to the band. It should 
be kept in mind that the 'coherent part' of the electron operator on the left side of Eq. (1) is not the full electron 
creation operator, as it does not contain terms that give rise to excited states of the boson degree of freedom [ p^ . 

It will be more useful to use hole operators rather than electron operators throughout this paper; we stress however 
that the discussion can be consistently carried out in electron as well as in hole representation. In terms of hole 
operators, the coherent part of the hole creation operator is 

4a ^[S+{T- S)n^.-,]c;J = S[l + Tn,^-a]c;J (2) 

Eq. (2) expresses the fact that the hole quasiparticle weight will increase with the local hole concentration, from S 
in the regime of low hole concentration to T for high hole concentration. The high degree of incoherence observed 
in high Tc cuprates for low hole doping implies << 1, and the fact that coherence is achieved for relatively small 
values of hole doping implies that the 'undressing parameter' 

T = I - 1 (3) 

is very large. T is the parameter that drives the transition to the superconducting state. Note that a large T 
necessarily implies S << 1, due to the constraint T < 1. For the normal state, Eq. (2) implies for the hole operator 

4a = s[i + ^r]c~J (4) 

with n the hole concentration per site, and for the hole number operator 

n,, = S^[l + ^T]^n^a = Z{n)h,„ (5) 

with Z{n) the hole quasiparticle weight. Eq. (5) implies that hole quasiparticles in the normal state become more 
coherent as the hole concentration increases. In the limit S* ^ 0, quasiparticles become completely incoherent in 
the normal state for low hole concentration and Fermi liquid theory breaks down. That limit is also described by 
the theory; in that limit, the transition to the superconducting state is a superconductor-insulator transition |p^ , |l6| . 
Even though for that particular limiting situation Fermi liquid theory does not describe the normal state, we stress 
that our approach is not a 'non-Fermi-liquid' approach, but instead is deeply rooted in Fermi liquid theory. 
Consider the bare kinetic energy in a tight binding model in terms of hole operators 

Hk^n = -J2*l^4aCj.+h.C.) (6) 

Replacement of the bare hole operators by the quasiparticle operators using eq. (2) yields 

Hhu = - XI (4aC]a + h.c.) (7a) 



= t^S^ll + T(n,,_, -I- nj- -t- T2n,,_,n,- (7b) 

Eq. (7) expresses the fact that the hopping amplitude of a hole quasiparticle will be increased, and as a consequence 
its kinetic energy will be lowered, as the local hole concentration increases; it is a direct consequence of the fact that 
the quasiparticle coherence increases with local hole concentration as described by Eq. (2) . For low hole concentration 
we can ignore the last term in Eq. (7b) and obtain 
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Hkin = ~ ^ [Uj + Atij{hi^^„ +nj^^cr)]{cla^jcr + h-c-) (8a) 



<IJ>,(T> 



Uj = S% (8b) 



Aty = Tty (8c) 



The kinetic energy of the form Eq. (8) is used in the model of hole superconductivity, and leads to pairing and 
superconductivity for low hole concentration in the presence of appreciable on-site and nearest neighbor Coulomb 
repulsion |9| . The condition for superconductivity to occur is 



T > v/(l + u)(l + w) - 1 (9) 

where u and w are dimensionless on-site and nearest neighbor Coulomb repulsions . Hence within this class of models 
superconductivity is intimately tied to increased quasiparticle coherence. Note that in a model with anisotropy Eq. 
(8) still implies 

^ = T (10) 

tij 

independent of direction. This assumption was used in our studies with the model of hole superconductivity and 
can be seen to be a necessary consequence of the fact that the At term in the Hamiltonian arises from quasiparticle 
undressing . A necessary consequence of Eq. (10) is that the superconducting energy gap function has the form 

Afc = A(efc) (11) 

and hence is constant over the Fermi surface (e^ = ep), even for an anisotropic band structure. Thus, Eq. (10) can 
be understood as a direct consequence of the undressing physics. Finally, Eq. (8) leads to superconductivity through 
kinetic energy lowering pO[ |. Hence, within the undressing scenario considered here, kinetic energy lowering as the 
system goes superconducting is intimately tied to s — wave symmetry of the superconducting order parameter as 
described by Eq. (11). 

III. GREEN'S FUNCTION: COHERENT PART 

The single particle Green's function is given by a sum of coherent and incoherent parts 

Gr.ir) = < Tcr^{r)cl^{0) >^ GZ\t) + GlT\r) (12) 

with T the time ordering operator. The coherent part of the Green's function is obtained by replacing the bare 
fermion operators in Eq. (12) by its coherent parts, given by Eq. (2) in terms of the quasiparticle operators: 

Glfir) = -S' < T[l + Tn.,i(r)]r:,t(r)][l + Tn^.^Wl^ltW] > (13) 
The normal and anomalous Green's functions for the quasiparticle operators 

Grsir) = - <Tcr^{T)dl^{0)> (14a) 



are given by the usual form ml 



FrsiT) = - <TCr^{T)£,l{0)> (14b) 



G(fc, zc.„) = . 4^ + . (15a) 



F{k,iuJn) = -UkVk[-. ^— '^—pT-] (15b) 

ILUn - Fk lU)n + J^k 
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where the coherence factors uu, Vk and quasiparticle energies are given by the usual BCS expressions 

= + (16a) 

-?. = ^(1-^) (16b) 

UkVk = ^ (16c) 



Ek = ^{ek-^^r + Ai (16d) 

and the gap function Ak is obtained from the BCS solution of the model of hole superconductivity 0, i.e. the kinetic 
energy Eq. (8) supplemented with on-site and nearest neighbor Coulomb repulsion. The single particle energy 
in these equations is given by Ck — Z{n)e1 5^(1 + nT)e^, with the bare kinetic energy given by the Fourier 
transform of {—t^j). 

It can be seen that the extra density operators in Eq. (13) will modify the normal Green's function introducing 
anomalous terms, similar to the anomalous terms that occur when calculating the expectation value of the kinetic 
energy Eq. (7) that lead to the optical sum rule violation [0. We expand Eq. (13), keeping only linear terms in the 
density as appropriate to the low hole concentration regime, and use mean field decoupling for the averages to obtain 

G™''(fc, zc^„) = S^[(l + nT)G{k, iLo^) + 2hTF{k, ilo^)] (17) 

with /o Ci|Ci| > the on-site pair amplitude in the superconducting state. We have also performed a space and time 
Fourier transform. It can be seen that the normal Green's function has acquired a contribution from the anomalous 
Green's function due to the density-dependent dressing. 

However the quasiparticle spectral weights derived from Eq. (17) are not positive definite and in fact can become 
negative in extreme parameter regimes. To remedy this we need to include higher order terms obtained from Eq. (13) 
by keeping terms with 6 fermion operators. Performing a similar mean field decoupling for these we finally obtain for 
the Green's function 

G™''(fc, za;„) = . ^\ + (1^^) 

lUJn — Ek lUJn + Ek 

with 

Zh = 5'[[1 + nT]uk - h^Vkf (18b) 
Z, = S''[[l + nT]vk + hTukY (18c) 

and 

/o =< 5aS.T ^ E ^(1 - 2/(i?fe))- (18d) 

The quasiparticle weights Zh and are clearly positive definite. Their sum is not conserved as function of density 
or temperature because of contributions from the incoherent part of the Green's function not contained in Eq. (13). 

In the absence of undressing (T = 0) the coherent Green's function Eq. (18) reduces to the usual BCS form 
except for the overall factor S"^ . In the presence of undressing (T > 0) Eq. (18) shows that the coherent part of 
the Green's function and spectral function will increase with hole density n, both for positive and negative energies. 
Furthermore, as the system goes superconducting the on-site pair amplitude /o develops a positive expectation value. 
From Eq. (18) this implies that the coherent spectral weight will decrease for positive energies (hole injection) and 
increase for negative energies (hole extraction) . The effect in the superconducting state will be largest for parameters 
where the on-site pair amplitude is large, which corresponds to short coherence length achieved in the strong coupling 
underdoped regime . The magnitude of these effects both in the normal and superconducting state depend on the 
magnitude of the undressing parameter T. We discuss the implications of these results in subsequent sections. 
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IV. RESULTS FOR QUASIPARTICLE WEIGHTS 



To illustrate the behavior emerging from the results of the previous section we consider now a specific example. 
The quasiparticle Hamiltonian is given by the kinetic energy Eq. (8) supplemented by on-site and nearest neighbor 
Coulomb repulsion 

Hcoui = U'^ fii^fiii + V fiiUj (19) 

i <ij> 

The BCS solution of this Hamiltonian 0] yields the quasiparticle energies 



Ek = ^{ek-fiy + Al = ^ a^ek - fi - + (20a) 

A,^A(efc) = A™(--^ + c) (20b) 

Ao = -A{p) (20c) 
a 

. = i|^A„ (20d) 

1 



(20e) 



' D/2' 

with A,„ and c parameters that depend on temperature and doping. The bandwidth D in these equations is given by 

D = Dh{l + nT) (21) 

with Dfi the bandwidth in the limit of zero hole concentration. The quasiparticle gap, i.e. the minimum quasiparticle 
excitation energy, is given by 

Eg = Ao (22) 

and occurs at momentum defined by 

4°'=M + ^ (23) 

However, if ej°' is below the bottom of the band, which occurs when the chemical potential is sufficiently below the 
bottom of the band at low hole concentration, Eq. (22) is not valid, and instead 



i^9 = y(-f -M)^ + A(-f) (24) 

We consider a two-dimensional square lattice with only nearest neighbor hopping and tij — th in Eq. (8). The 
quasiparticle bandwidth as the hole concentration goes to zero is Dh — "izth, with z — A the number of nearest 
neighbors to a site. We choose parameters 

D,, = 0.2eV 
U = beV 
V = OMeV 

T = 19.2 (25) 

which imply At = TDh/2z = OASeV. For the present purposes we need not specify the magnitude of the parameter 
5^, which determines the relative weight of coherent and incoherent contributions to the spectral function. 
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These parameters yield a maximum Tc versus hole concentration of T^"'^ = 94K, as shown in Fig. 1(a), for optimal 
doping n ~ 0.045. The minimum quasiparticle excitation energy at low temperatures is shown in Fig. 1(b). At low 
hole concentration it does not go to zero as Tc does because the chemical potential falls below the bottom of the band 
and Eg is determined by Eq. (24) rather than by Eq. (22). The behavior of the chemical potential and the bottom 
of the band versus hole concentration is shown in Fig. 2. The chemical potential crosses the bottom of the band at 
n ~ 0.038, and e° (Eq. (23)) crosses the bottom of the band at n 0.034. 

The on-site pair amplitude /o that enters in the expressions for the quasiparticle weights is shown in Fig. 3. As 
function of doping it follows approximately the behavior of the critical temperature and of the gap parameter Aq (not 
shown). At low hole doping it goes to zero because the carrier concentration goes to zero, at high hole doping it goes 
to zero because the coherence length is diverging [ p^ . As function of temperature, /o behaves approximately like the 
gap, going to zero at as (T^ — T)^/^. 

Next we consider the behavior of the quasiparticle weights Ze and as function of temperature. Figure 4 shows 
the results at the (normal state) Fermi energy, — fi, for the optimally doped case (n = 0.045). The values are 
normalized so that Z^ and Zh would be 0.5 for T = 0. The dashed line shows the value the weights would have 
for /o — 0: it is temperature independent and larger than 0.5 because of the undressing due to the average carrier 
concentration n. The effect of onset of superconductivity is to increase Z^ as the temperature is lowered and to 
decrease Z^,. This indicates that there is extra amplitude for electron creation, and less amplitude for hole creation. 
It may thus be interpreted as a shift of the chemical potential as superconductivity sets in, giving increased hole 
occupation as the temperature is lowered, or equivalently a shrinking of the electron Fermi sea. This is a surprising 
result of this calculation, and its implications will be discussed in subsequent sections. Note that the weight Ze 
increases by almost a factor of 2 between T = Tc and T = 0. The magnitude of the increase of course depends on the 
magnitude of the undressing parameter T , and would be larger or smaller for larger or smaller values of T respectively. 
By adjusting the values of on-site and nearest neighbor Coulomb repulsion in the model it would be possible to obtain 
the same maximum Tc with different values of T, as discussed in previous work Nevertheless we believe that the 
parameters chosen for this example may be representative of the situation in high Tc materials. 

Note also that the total weight of the spectral function Ztot = Zc + Zh increases as the temperature is lowered 
below Tc- This indicates that overall there is more coherence in the superconducting than in the normal state, in 
accordance with qualitative expectations, and this extra spectral weight is transfered from the high energy incoherent 
part of the spectral function as will be discussed in the following section. However part of the enhancement of Zc at 
low temperatures relative to its value at Tc can be attributed to spectral weight being transfered from negative to 
positive energies (i.e. corresponding depletion of Z^) in addition to spectral weight transfer from the incoherent part 
of the spectral function. 

Similarly Fig. 5 shows the results for an overdoped case, n = 0.1, with Tc = 68K. The behavior is qualitatively 
similar as the optimally doped case, however the effect of the onset of superconductivity on the spectral weights is 
considerably smaller because the system is already more coherent in the normal state. This is indicated by the larger 
values of all the spectral weights relative to the values of the case shown in Fig. 4, due to the enhanced coherence 
arising from the increased hole concentration. For a much higher hole concentration, as Tc approaches zero, the 'gap' 
between Zc and Zh in the superconducting state closes, as shown in Fig. 6. It remains always nonzero however as 
long as Tc is nonzero, and there is always some spectral weight transfer from the incoherent region as long as Tc is 
nonzero. 

Next we consider the spectral weights for other values of momentum. Figure 7 shows results for — /i > 0. Recall 
that we are using hole representation, so efe — /i > means inside the filled Fermi sea for electrons. In the normal 
state Ze — 0: since the electron state is full, no new electron can be created in it. Just as in the conventional 
BCS case, as superconductivity sets in the state becomes partially occupied and Ze ^ , and Zh correspondingly 
decreases. However, unlike conventional BCS, Ze and Zh cross in our case and at low temperatures the weight for 
creating an electron is larger than that for creating a hole, even though we are inside the filled normal state Fermi 
sea. Clearly this implies that the chemical potential in the superconducting state has changed. Figure 7(b) shows 
that for Ck — ^l = IQ.lmeV the weights for electrons and holes coincide at low temperatures; this momentum then 
corresponds to the new Fermi momentum, k'p, in the superconducting state. For even larger efc, Zc becomes smaller 
than Zh as in the conventional case, as shown in Fig. 7(c). 

The behavior for negative energy (outside of the electron Fermi sea) is shown in Fig. 8. Here, Cfc — was chosen 
to be at the bottom of the hole band in the optimally doped case. The weight for electron creation is much larger 
than in the conventional case. Note also that Ze first decreases and then increases as the temperature is lowered, 
and as r — > it becomes even larger than its value in the normal state. Such a situation, which is never seen in the 
conventional case, is possible here due to the non-conservation of Ztot because of the transfer of spectral weight from 
the incoherent part to the coherent part of the spectral function as the temperature is lowered. 

Figure 9 shows the spectral weights for an overdoped case, n = 0.1, for values of the momentum above the electron 
Fermi surface, at the Fermi surface and below the Fermi surface. The behavior is qualitatively similar to that for 
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the optimally doped case, although the differences between the conventional and our case are less pronounced here 
because this is a weaker coupling regime. In figure 10 we show the spectral weight for an underdoped case, n = 0.02. 
Here the chemical potential is below the bottom of the band, so the situation comparable to Fig. 4 cannot be attained. 
Fig. 10 shows the behavior of the spectral weight for at its lowest possible value, the bottom of the band, which is 
qualitatively similar to other cases where tk is above jj. such as Fig. 7(a). 

Next we consider the behavior of the quasiparticle weights at the chemical potential versus doping in Fig. 11. The 
upper dot-dashed line is the total spectral weight in the superconducting state, and the dotted line below it is the total 
spectral weight in the normal state. The difference between the two is the spectral weight transfered from high energy 
incoherent processes as the system goes superconducting; this difference approaches zero in the overdoped regime. 
The full lines denote the quasiparticle weights in our case, the dashed lines the usual BCS results (u^ = Zh,v1 = Z^), 
which increase approximately linearly with n due to the normal state increased coherence with doping and are equal 
for efc = fi. For low dopings however the chemical potential falls below the bottom of the band and hence we take 
Efc at the bottom of the band rather than at fi, this is why the two dashed lines diverge at low dopings. In our 
case, the weight for electron creation (solid curve labeled Z^) is seen to increase rapidly with doping and then taper 
off for high doping; this latter effect is due to the reduction of the on-site pair amplitude /o for high doping as the 
coherence length becomes large . The quasiparticle weights Z^ and Zh approach each other and the BCS value for 
high doping, as expected. Note also that there is a narrow doping regime where the electron weight Zg is even larger 
than the total weight in the normal state (dotted line). This situation can never occur in the conventional BCS case. 

We believe the behavior exhibited by Z^, in Fig. 11 is relevant to the understanding of the angle resolved photoemis- 
sion results discussed by Ding et al ||l| . In their work, the quasiparticle weight in the superconducting state extracted 
from photoemission spectra shows similar qualitative behavior to the behavior exhibited by Z^^ in Fig. 11. We will 
discuss the relation between Z^ and the experimental quantity in a subsequent section. Ding et al also plot ZA^ the 
product of their extracted quasiparticle weight and the gap inferred from the photoemission spectra, and point out 
that its behavior rougly follows the bell-shaped curve of Tc. Our calculation shows similar behavior, as shown in Fig . 
12. Note that the quasiparticle gap itself remains finite in our calculation as the hole concentration goes to zero |l^] 
(Fig. 1 (b)), and also experimentally 



V. GREEN'S FUNCTION: INCOHERENT PART 



To calculate the incoherent part of the Green's function we now consider a specific model, the generalized Holstein 
model discussed in I. Our calculation follows closely the calculation of Alexandrov and Ranninger [|l9| for the conven- 
tional Holstein model, and we refer the reader to their seminal work for details which are common to both situations. 
The site Hamiltonian for our case is given by [p^ 

H — huJaa^a + gfiujQ{a) + a)(n| + — ^n'\ni) -\- Un^ni (26) 

The case of Alexandrov and Ranninger corresponds to 7 = 0. Using a generalized Lang-Firsov transformation [T^ ] 
the quasiparticle (polaron) operators cia are related to the bare fermion (hole) operators by 

C-icr — CicrXicr (27a) 



X,^ = g-3(at-a.)(l-7ft.-.) (27b) 

In contrast to Eq. (2), the operator da here is the full hole destruction operator, including coherent and incoherent 
parts. The coherent part results from the expectation value of the X-operators in the zero boson subspace, 

< >= e-4(i-T».-) (28) 

and in particular 

X,^{n,^^ =0) = e-"^ ^ S (29a) 



T 



(29b) 



in accordance with Eq. (2). 
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We wish to calculate the Green's function 
G(m,r) = -<rcoT(r)4^(0) 

— „ < yg9(aJ(r)-ao(r))(l-7fio4(T))~^^|-^-j~t ^ j-Q-jg-ff(at^ (0) -a™ (0)) (l-7fi„i (0)) ^ ^gg-j 

We expand the exponentials in Eq. (30) using the operator relation 

and decouple averages over bosons and fcrmions, following Alexandrov and Ranninger. This leads to 

r) = <JQ{m,T) < -Tco|(r)c|„-^(0) > 

+ [(cri("i,T) -ao{m,T)][< -Tnoi(r)coT(T)c|„^(0) > + < -rcoT(T)c|'„^(0)n„j.(0) >] 

+ [CT2(m, t) - 2(71 (m,r) + (To(m, r)] < -Tnoi(r)co|(T)cJ„^(0)n™|(0) > (32) 

The boson Green's functions are defined as 

a,+j(m, r) =< e9(l-7)*(a;(r)-ao(r))g-5(l-7)^(aL(0)-a,„(0)) > ^33^ 

with i, j = 0, 1. At low temperatures they are given by 

i„(m, t) = S'-'T'll - 4„.,„ + «.,..oe-»'<'-''"'"''] (34a) 



and in frequency space by 



(m, *c.„) = 52-"T" [,5„,o/3 + <5,„,o V „ ] (35) 

We next decouple the fermion averages with the same mean field procedure used to calculate the coherent part of the 
Green's function , and calculate the Fourier transform 

G(fc, iun) = f dre'^"^ V e''"'Go,n{r) (36) 

The complete Green's function is 

G(fc, iUn) = Gcohik, iuJn) + Ginc{k, iUJn) (37) 

For each term of the coherent Green's function, Eq. (18), there is a corresponding term in the incoherent Green's 
function. All terms in the coherent Green's function are of the form 

G:„-(fc,^^„) = c5^-"r° (38) 
with < a < 2 and s = +/ — 1. The corresponding term in the incoherent Green's function is 

'""^ N ^ li^n - s{Ek' ~ IlUo) ILUn - s{Ek' +IuJq)' 

The spectral function is obtained as usual from 

A{k,LL!) = -ImG{k,iu!n ^ to + iS) (40) 

and results from Eqs. (18), (37)-(39). In particular, the lowest order normal part of the incoherent spectral function 
is given by 



A^^,{k, a;) = 52 X 5^ ^[1 + n[(l - 7)'e^^ - 1]] 



fc' 

,,2 



+ i;^, K.(5(w - lu>o + Ek') + {1- nk')S{uj + Iujq + Ew)]] (41) 
and the lowest order anomalous contribution by 



A1^,{k, a;) = 52 X 2/0 X 5^ ^ [1 + n[(l - 7)'e^« - 1]] 

1=1 

X — ^^{-Uk'Vk') X [(1 - nk')6{uj - lujo - Ek') + nk'5{uj + Iloq - Ek') 

k' 

- nk'S{w - lujo + Ek') - (1 - nk')S{w + Icoq + Ek')]. (42) 



VI. RESULTS FOR THE SPECTRAL FUNCTION 

The spectral function for the models considered here is of the form 

A{k, uj) = Acoh{k, w) + Ainc{k, w) (43a) 

Acoh{k, oj) = Zh5{oj - Ek) + ZeS{u + Ek) (43b) 

A^ncik, uj) = -ImGincik, UJ + iS) (43c) 

where the quasiparticle weights Zh and Z^,, and the incoherent Green's function Ginc, were discussed in the previous 
two sections. As seen in Sect. IV, the quasiparticle weight acquires a positive contribution from the onset 
of superconductivity. In a spectroscopic experiment usually only one side of the spectral function is sampled, as 
the other side is suppressed by the Fermi function. The quantity that will display the enhanced coherence due to 
undressing exhibited by Z^ is 

Io{k,ij) = A{k,uj)f{uj) (44) 

with / the Fermi function. In an experiment there will typically be broadening from experimental resolution, which 
results in 

I{k,u) = j duo'F{uj -uj')Io{k,u') (45) 

being measured, with Fiuo) a Gaussian with width a^^. There could also be other sources of broadening of the 5- 
functions in the expressions Eq. (43) from lifetime effects. Just as the spectral function, the measured spectrum Eq. 
(45) will have coherent and incoherent contributions 

7(fc, w) = Icoh{k, w) + Iinc{k, w) (46) 

arising from the coherent and incoherent parts of the spectral function respectively. 

Figure 13 shows results for the coherent spectra at the Fermi energy for an underdoped (n = 0.02, labeled ud), 
optimally doped (n = 0.045, labeled op) and overdoped (n = 0.1, labeled od) case for the parameter values used in 
Sect. V. For the underdoped case with the chemical potential below the bottom of the band the value of ek at the 
bottom of the band was used. The dashed lines show the spectra in the normal state at Tc, and the full lines in 
the superconducting state at T = O.lTg. For each doping, as superconductivity onsets the peaks shift to the left due 
to the opening of the superconducting gap. Furthermore, the peaks grow in magnitude due to the behavior of Z^, 
discussed in Sect. IV. As function of doping the peaks grow in magnitude both in the normal and superconducting 
state due to the enhanced coherence with increased number of carriers. The superconducting peak in the od case is 
shifted to the right with respect to the op case because the superconducting gap is smaller in that region (Fig. 1(b)). 
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When we include the incoherent part of the spectra the smaller normal state peak can become almost invisible. 
The results will depend of course on the specific parameters chosen to describe the incoherent background, and we 
are not suggesting that we are in a position to determine them from first principles. In Figure 14 we show results for 
a particular set of parameters for the generalized Holstein model. In addition to the parameters already discussed in 
Sect. IV, including the value of T, the new parameters needed are S^, luq and a broadening factor, given in the figure 
caption. Note that in the underdoped case (a) the peak in the normal state has become almost invisible, while a 
sharp peak and a dip are seen in the superconducting spectrum. The dip arises because the background term arising 
from the second term in Eq. (39) for ak> = u^., , s = — 1, is pushed to more negative energies as the superconducting 
gap opens. As the doping increases the normal state peak becomes more visible, and the overdoped case shows more 
conventional behavior. Note that the scale in the figures changes with doping and the magnitude of the peaks increases 
with doping. 

Figure 15 shows the temperature dependence of the spectra for the overdoped case. The normal state peak is 
pushed back continuously as the superconducting gap opens up. In addition, for our case (a) the peak grows in 
magnitude. To highlight the difference with conventional BCS theory we show in Fig. 15(b) what is obtained with 
the same parameters in the absence of the term /o in Eq. (18). The peak here first becomes lower and then increases 
again as the temperature is lowered, but it is always lower or equal to the normal state peak. It is easily seen from 
the BCS formula that this property is generally true also for other values of the momenta. 

Similarly figure 16 shows the temperature dependence of the spectra in the underdoped case. Here, rather than 
the peak moving continuously, a new peak grows in the superconducting state. The presence of two peaks has not 
been seen experimentally in photoemission to our knowledge, possibly because of experimental resolution. For the 
BCS case (b) the peak in the superconducting state is much smaller than in the normal state, while for our case (a) 
the opposite is true. Results for the temperature dependence in the optimally doped case show behavior intermediate 
between the overdoped and underdoped cases shown. 

VII. THE CUPRATES 

We have seen in the previous sections that in systems where superconductivity arises from undressing there is a 
signature of the formation of the condensate in the single particle spectral function. Specifically, it arises in Eq. (17) 
from the term involving /q, the on-site pair amplitude. Ding et al ^ and Feng et al discussing experimental results 
of angle resolved photoemission in cuprates have recently emphasized precisely that feature of the observed spectra, 
and correlated the growth of the peak in photoemission to quantities related to the superconducting condensate 
such as the superfluid density and the condensation energy. The spectra calculated within the present theory in the 
previous section resemble in several aspects the experimental observations in photoemission along the (tt, 0) direction. 

Unfortunately, as the alert reader has undoubtedly noticed, the results presented in the previous section with nega- 
tive a; in a hole representation correspond to hole destruction, or electron creation, that is, inverse photoemission. 
It is for that case that the experiment would sample Z^,, the quasiparticle weight for electron creation. Instead, if 
we calculate spectra for direct photoemission we would find that quasiparticle peaks are suppressed by the onset of 
superconductivity due to the behavior of Zh discussed in Sect. IV. 

The present theory does not allow for a switch in the role of the weights Z^ and Zh'- electron- hole asymmetry, of 
the sign assumed here, is central to the theory. Does this then imply that the theory is irrelevant for description of 
the cuprates? 

We believe this is not the case. We propose that in fact, the photoemission experiments along the (tt, 0) direction 
close to the (7r/a,0) point sample the part of the spectral function discussed in the previous section, corresponding 
to hole destruction, or electron creation. 

How can photoemission sample electron creation? Recall that in the theory of hole superconductivity the relevant 
orbitals are oxygen pn orbitals in the planes There are however also the oxygen pa orbitals, strongly hybridized 
with the Cu d^^-y^ orbitals. Suppose that in a photoemission experiment along the (tt, 0) direction the largest matrix 
element couples to destruction of a d^2_y2 electron. This would not directly couple to the band responsible for 
superconductivity; however, that process could induce the destruction of an oxygen hole in the pn orbitals. The 
proposed situation is schematically depicted in Fig. 17, in an electron representation. Before the photon comes in 
there is one electron in each Cu'^^ atom neighbor to a given O atom, and two holes in the pn orbital on that O atom. 
We assume the energy level structure shown in Fig. 17: an electron from Cu^+ cannot 'fall' onto the O pn orbital 
because it is Coulomb-repelled by the electron in the other Cu atom. When a photon comes in and knocks out one 
of the electrons in a Cu, the other electron can 'fall' onto the O pir orbital, thus destroying an O hole and sampling 
the quasiparticle weight Z^. 

It is clear that this qualitative explanation needs further elaboration and experimental confirmation to be convincing. 
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Nevertheless we also point out that it suggests an explanation for why the sharp peaks seen in photoemission along 
the (tt, 0) direction are not seen along the (tt, tt) direction ||2l|] : since Cu dx2-y2 orbitals point along the principal axis 
in the planar square lattice, the coupling to the photon along the (tt, tt) direction is likely to be much smaller. For that 
direction the larger coupling may be to the O pn band itself, in which case inverse rather than direct photoemission 
would show the enhanced coherence. It is possible that some indication of this effect may have already been seen in 
tunneling experiments [^2| . 



VIII. CONCLUSIONS 



In this paper we have continued exploring the consequences of the physical principle proposed in I: that, in at least 
some electronic materials in nature, the dressing of quasiparticle carriers is a function of the local carrier concentration, 
and becomes smaller as the local carrier concentration increases. This physical principle leads to superconductivity 
occuring in these systems because of lowering of the carriers kinetic energy upon pairing. The superconducting 
transition, and many of the features of the superconducting state, have already been discussed earlier within the 
theory of hole superconductivity . 

In this paper we explored the consequences of this principle for the single particle Green's function in the supercon- 
ducting state. The central result of this paper, Eq. (18), demonstrates that formation of the superfluid condensate 
will influence the behavior of the single particle spectral function. Eq. (18) is thus the generalization of the BCS 
spectral function for systems where superconductivity is driven by undressing. Surprisingly, the results show that the 
enhanced coherence in the superconducting state is displayed in the quasiparticle weight for electron creation but not 
for electron destruction. We calculated the behavior of the quasiparticle weights as function of temperature, doping 
dependence and momentum, and highlighted their differences with conventional BCS theory. 

Furthermore, we discussed the calculation of the full spectral function including the incoherent contribution for one 
particular model where superconductivity occurs through undressing, a generalized Holstein model. Our calculation 
was performed within the Lang-Firsov approximation |19[| , and it should be interesting to see whether the qualitative 
results survive a more exact treatment |2^j2^ . 

Results for the full spectral function showed several features that resemble experimental observations in photoemis- 
sion experiments in high Tc cuprates ||l|,^, in particular enhanced coherence, as displayed by the quasiparticle peak 
in the spectra, when the system enters the superconducting state and as the carrier concentration increases both in 
the normal and in the superconducting state. 

This study was strongly motivated by the beautiful experimental results and insightful analysis of the photoemission 
experiments [^]j2j. Thus it is perhaps disappointing that at the end of the day our calculation predicts these effects, 
in the simplest one-band model, arising in inverse rather than in direct photoemission. Thus some readers may 
conclude that our calculation is not more than an academic exercise. However, as discussed in Sect. VII, we believe 
there is a plausible scenario by which the spectral weight for electron creation would be sampled in the photoemission 
experiments in the cuprates. 

While the theory discussed here predicts s-wave rather than d-wave superconductivity we believe it is remarkable 
how many of the features that appear to be part of the phenomenology of high Tc cuprates it exhibits, as a consequence 
of the single assumption of a large value of the undressing parameter T: (1) incoherence in the normal state at low 
hole concentration; (2) increased coherence with doping in the normal state; (3) transition to superconductivity for 
low doping, dissappearing for high doping, and bell-shaped Tc versus hole concentration; (4) increased coherence 
as the system goes superconducting; (5) superconducting transition driven by kinetic energy lowering, optical sum 
rule violation; (6) non-decrease of the quasiparticle gap at low hole density when Tc is going to zero. This latter 
feature arises in our model from the fact that as the hole concentration decreases and the band becomes narrower 
the chemical potential falls below the bottom of the band ||l^ ; we believe that many of the unusual properties of 
underdoped cuprates follow from this simple fact, and in particular that the observed pseudogap is simply the energy 
difference between the bottom of the band and the chemical potential. [ p5p^ ] 

If the theory of hole undressing discussed here describes the cuprates it is likely that it is more generally applicable, 
because it is based on very general principles. In this regard we note that one of the paradoxes of the conventional 
explanation of superconductivity is that it is thought to originate in an electron-boson (the electron phonon) coupling 
that opposes conductivity, i.e. gives rise to resistivity, in the normal state. In a sense the present theory eliminates 
this paradox. Coupling to a boson is certainly necessary, and that coupling gives rise to enhanced resistivity in the 
normal state due to enhanced effective mass, but superconductivity arises from a process whereby the coupling to 
that boson is reduced as carriers pair and the system becomes superconducting. The old paradox is however replaced 
by a new one, that in order for heavily dressed 'confined' carriers to become less dressed, or 'freer', it is necessary for 
them to bind in Cooper pairs. 
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We also note that the principle on which the present theory is based, that an increase in the local hole occupation 
causes undressing, is likely to be more general than as expressed by Eq. (2): rather than just by the same site 
occupation, undressing may also be enhanced by hole occupation of neighboring sites, and also by neighboring bond 
occupation. The possible implications of this for superconductivity and other instabilities of metals will be discussed 
in future work. 

If indeed the essential physics of high Tc is hole undressing, what makes a material a high Tc superconductor? 
Presumably, the fact that quasiparticles are heavily dressed in the normal state together with the fact that the 
undressing process that occurs when the local carrier concentration increases is particularly efficient. Both of those 
facts are necessary conditions for high superconductivity by giving rise to a large T parameter. We will not discuss 
here what aspects of the chemistry of the cuprates would favor this situation p7[ . However, conversely, we may 
conclude that the reason for a material not being a high superconductor would be a small value of the parameter 
T, either because quasiparticles are not heavily dressed in the normal state (e.g. the case of Aluminum), or, because 
the quasiparticle dressing in the normal state may not be strongly dependent on the local carrier concentration (e.g. 
the case of 'heavy fermion' systems). 
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FIG. 1. Superconducting transition temperature versus hole doping n, number of holes per planar oxygen, for parameters 
given by Eq. (25). (b) Minimum quasiparticle excitation energy at low temperatures versus doping. 

FIG. 2. Chemical potential /x at low temperatures, and band bottom (— D/2) versus hole doping for the parameters Eq. 
(25). The chemical potential falls below the bottom of the band for hole concentration n ~ 0.038. For fixed hole concentration 
fj, increases as the temperature is lowered above Tc, particularly for low hole concentration, and stays approximately constant 
below Tc for all hole concentrations. 



FIG. 3. On-site pair amplitude /o (a) versus hole doping at low temperatures, and (b) versus temperature at optimal doping. 

FIG. 4. Quasiparticle weights at the Fermi energy (e^ — fi) versus temperature for the optimally doped case. Tc = 94_R". 
Ztot is the sum of Zh and Ze. The corresponding BCS results are equal to each other and independent of temperature (dashed 
line). 



FIG. 5. Same as Fig. 4 for an overdoped case, n = 0.1, with Tc = 68K. 



FIG. 6. Same as Fig. 4 for a highly overdoped case, n = 0.2, with Tc = 2AK. 

FIG. 7. Spectral weights at optimal doping versus temperature for momentum inside the electron Fermi surface. The dashed 
Unes give the BCS values (Eq. (18) with /o = 0, ~ Zh, v't ~ Zh). The upper dashed hue corresponds to ul, the lower one 
to v^. 

FIG. 8. Spectral weights at optimal doping versus temperature for momentum outside the electron Fermi surface. The 
dashed lines give the BCS values, the upper dashed line corresponds to v^, the lower one to u^.. 



FIG. 9. Spectral weights versus temperature for an overdoped case, n = 0.1, and momentum (a) outside, (b) at, and (c) 
inside the electron Fermi surface. 

FIG. 10. Spectral weights versus temperature for an underdoped case, n = 0.1, and momentum inside the electron Fermi 
surface. The value of — /i = 6.5meV corresponds to at the bottom of the hole band. 

FIG. 11. Quasiparticle weights versus doping at low temperatures. The weights are computed at the chemical potential 
when it is inside the band, and at the lower hole band edge in the underdoped regime when the chemical potential is below 
the band edge. The lower dashed line gives the BCS values m| and v^, which are equal when fj, is inside the band and separate 
into two (upper corresponding to lower to vl) when /i is below the hole band edge. The dot-dashed line Ztot gives the total 
weight Zh + Ze and the dotted line close to it the corresponding BCS total weight. Note that Ze rises approximately linearly 
with doping for low hole doping and levels off in the overdoped regime. All results approach the BCS values for high doping as 
fo approaches zero, but remain different from the BCS values as long as Tc is nonzero. 

FIG. 12. Product of electron quasiparticle weight Ze and minimum quasiparticle excitation energy Eg at low temperatures 
versus hole doping. The dashed line gives Tc versus doping. Note that Ze x A peaks at somewhat higher values of hole doping 
than Tc does. 

FIG. 13. Coherent part of the spectral function multiplied by the hole Fermi function /{lo) and broadened by a Gaussian 
function, Eq. (43), with au, = 5meV. The dashed lines give the results at Tc, the full lines the results at low temperatures 
(T = O.lTc). od, op and ltd denote overdoped (n = 0.1), optimally doped (n = 0.045) and underdoped (n = 0.02) regimes. 
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FIG. 14. Results for the full spectra Eq. (43) at T = Tc (dashed lines) and at low temperatures, T = O.lTc (solid lines). 
The incoherent part of the spectrum was modeled with a generalized Holstein model with wo = 5meV and gaussian broadening 
for the 5 functions a — 15meV. The band narrowing parameter is S'^ = 1/500, and T = 19.2, corresponding to g — 2.49 and 
7 = 0.45 in Eq. (26). The momentum is given by the normal state Fermi momentum for the optimally doped and overdoped 
cases, and by the momentum corresponding to the bottom of the hole band for the underdoped case. 

FIG. 15. Spectra in the overdoped case n = 0.1 at the chemical potential for various temperatures. The dashed, dotted, 
dot-dashed, dotted and full lines, with the peak moving towards the left, correspond to T/Tc = 1, 0.9, 0.8, 0.6 and 0.1 respectively. 
In (b) the corresponding results for the BCS case, taking /o = in Eq. (18), are shown. Note that the peaks in the 
superconducting state are always lower than that in the normal state in the BCS case. 

FIG. 16. Same as Fig. (15) for an underdoped case, n = 0.02, for momentum corresponding to the bottom of the hole band. 
The same line convention as in Fig. (15) is used. 

FIG. 17. Schematic proposed explanation of how electron extraction in a photoemission experiment can give rise to hole 
destruction in the oxygen band responsible for superconductivity. The Cu energy level arises from a hybridized Cu d^2_y2-0 
pa orbital, the O energy level corresponds to an O pn planar orbital. In (a), an electron in the Cu d^2_yi orbital is knocked 
out by an incoming photon; the electron from a neighhboring d^2_y2 orbital then falls onto the O pi: orbital (b), destroying an 
O hole. 
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